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The Hilbert-Polya conjecture states that the imaginary parts of the zeros of the Riemann zeta 
function are eigenvalues of a quantum hamiltonian. If so, conjectures by Katz and Sarnak put this 
hamiltonian in Altland and Zirnbauer's universality class C. This implies that the system must have 
a nonclassical two- valued degree of freedom. In such a system, the dominant primitive periodic orbits 
contribute to the density of states with a phase factor of —1. This resolves a previously mysterious 
sign problem with the oscillatory contributions to the density of the Riemann zeros. 



o 

(N 

> 
O 



43 ■ 

43 



> 
(N 

in 
o 



The distribution of prime numbers among the integers 
is a fundamental problem of number theory (see e.g. 
It is closely connected to the properties of Dirichlet L 
functions (including the Riemann zeta function) , defined 



L(s,x) ■-- 
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for Re s > 1 and by analytic continuation elsewhere, 
where x( n ) is a primitive Dirichlet character: x( n ) is 
periodic with smallest period d, has magnitude one or 
zero, is zero if d and n are not coprime, and obeys 
x(mn) = x( m )x( n )- The Riemann zeta function is given 
by C( s ) := L(s,l). According to the generalized Rie- 
mann hypothesis, any zero of L(s,x) with < Res < 1 
is on the critical line Re s — | ; these are the nontrivial 
zeros, which we will write as pk = \ + ijk- The gener- 
alized Riemann hypothesis implies that each 7^ is real, 
and this in turn can be shown to imply that the num- 
ber of primes less than x in the arithmetic progression 
a,a + d,a + 2d, . . . (with a less than and coprime to d) is, 
in the limit of large x, 
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for all e > 0, where ip(d) is the number of integers less 
than and coprime to d (the Euler totient function), and 
hi(x) is the logarithmic integral function. The exponent 
of x in the error term increases to ^ + max Im 7^ + e if 
the generalized Riemann hypothesis is false. 

It is an old idea, now generally known as the Hilbert- 
Polya conjecture (see Q for a historical review) , that the 
nontrivial zeros of each L function are the eigenvalues of 
an operator (on some Hilbert space) that takes the form 
i + iH , where H is self-adjoint; each L function would 
have a different H. Furthermore, L(i + iE,x) is con- 
jectured to be proportional to the spectral determinant 
det(E — H); since the eigenvalues of a hermitian oper- 
ator must be real, the Hilbert-Polya conjecture implies 
the generalized Riemann hypothesis. An explicit con- 
struction of the ff's for the different L functions (or, less 
ambitiously, just for the Riemann zeta function) would 
therefore be extremely important. 



A large body of analytic and numerical work strongly 
supports the Montgomery-Odlyzko law (see e.g. |3J), 
which states that the statistical distribution of the 7fc's 
for each L function is the same as the Wigner-Dyson 
distribution of the eigenvalues of large hermitian matri- 
ces with real diagonal entries and complex off-diagonal 
entries, each selected from a gaussian distribution; this 
is the gaussian unitary ensemble (GUE) A large 

body of analytic and numerical work also strongly sup- 
ports the Bohigas-Giannoni-Schmit conjecture [5(, which 
states that the energy eigenvalues of the hamiltonian for 
a system that is classically chaotic, and not time-reversal 
invariant, also obey the GUE distribution. This leads to 
the generalized Berry conjecture [fj|: the operator H for 
each L function can be obtained by quantizing a classi- 
cally chaotic system that is not time-reversal invariant. 

Katz and Sarnak Q have conjectured that L func- 
tions corresponding to Dirichlet characters that are real 
ix( n ) = 0,±1] and even [%(— 1) = +1] form a "family" 
(that includes the Riemann zeta function) whose mem- 
bers are related (in some fashion) by a symplectic sym- 
metry, and furthermore that the spacings of the 7fc's for 
each member of this family is governed by the distribu- 
tion of eigenphases of random unitary symplectic matri- 
ces. This agrees with the GUE distribution for 7/. 3> 1, 
and predicts a gap in the spectrum near zero; this is well 
supported by numerical evidence from these L functions 
0, H[ . Other proposed families have unitary or orthogo- 
nal symmetries. 

Altland and Zirnbauer [9( have classified the possible 
symmetry classes of quantum hamiltonians. The distri- 
bution of 7fc 's found by Katz and Sarnak is a predicted 
property of the energy eigenvalues for a chaotic system 
in Altland and Zirnbauer's class C. We therefore inter- 
pret the Katz-Sarnak conjecture, in the context of the 
Hilbert-Polya conjecture, to mean that the quantum sys- 
tem corresponding to the Riemann zeta function (or any 
other member of its symplectic family of L functions) 
should have a hamiltonian in class C. 

A hamiltonian in class C takes the form of a generator 
of USp(N); more specifically, 

H = A + a-S, (3) 
where A is a hermitian operator that (when expressed 
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as a matrix in a suitable basis) is imaginary and anti- 
symmetric, and each Si (i = 1,2,3) is a hermitian op- 
erator that (when expressed as a matrix in the same 
basis) is real and symmetric; finally, <Zj is a Pauli ma- 
trix acting in an additional two-dimensional Hilbert 
space. This extra "nonclassical two-valuedness" ( "klas- 
sisch nicht beschreibbare Zweideutigkeit" , Pauli's [lOj de- 
scription of electron spin) is a previously unrecognized 
essential ingredient in any attempt to construct a quan- 
tum hamiltonian with eigenvalues corresponding to the 
imaginary parts of the nontrivial Riemann zeros. 

Next, consider the "completed" zeta function A(s) := 
r 0O (s)C(s), where r oo (s) := 7i-- s / 2 r(s/2) and T(z) is the 
Euler gamma function. The completed zeta function 
obeys Riemann's functional equation A(s) = A(l — s), 
and is real on the critical line; the zeros of A(s) coin- 
cide with the nontrivial zeros of C( s )- It follows that the 
number of zeros of £(s) on the critical line with imaginary 
part between zero and E > is given by 



N(E) = i Imlog A(± +e + iE) + l, 



(4) 
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ere e is a positive infinitesimal (ill ] - We can write 
N(E) as the sum of a smooth contribution and an oscil- 
lating contribution [6|: 

N{E) =N(E) + N osc (E), (5) 
N(E) =ilmlo g r oo (i+i J B) + l 

= IMI)-| + l+« (6) 
JV OBC (B) = ilmlogC(| + e + iS). (7) 

Using the Euler product formula £(s) = Yl p (l — p _s ) _1 , 
where p is a prime, we get the formal expression 

N osc (E) = -ilmVlogOL -p-CVa+u*)) 
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This expression is formal because the Euler product does 
not converge on the critical line. Its value is in its sim- 
ilarity to the corresponding expression for the number 
of energy eigenvalues less than E of a hamiltonian for a 
classically chaotic system whose classical periodic orbits 
are all isolated and unstable. For a system without the 
two- valued quantum degree of freedom required by class 
C, the smooth contribution is given by the Weyl formula 
(see e.g. 0) 



N(E) 



d^x d?p 



9(0 < h(x,p) < E), 



(9) 



where 6(5) = 1 if S is true and if S is false, / is the 
number of classical degrees of freedom, and h(x,p) is the 



classical hamiltonian 13j. The oscillating contribution 
is given by a formal sum over primitive periodic orbits 
(labelled by po) and their repetitions (labelled by r) , 



po r 



-r|det(iW>-J)|i/2' 



where the primitive orbit has action S po (E), Maslov 
phase /j, po (E), and stability matrix M po (E). 

If we hypothesize a dynamical system in which the 
primitive periodic orbits are labelled by prime numbers 
[6j , then eq. (flQ|) bears a strong resemblance to eq. (jHJ ■ 
However, there are two well known problems with get- 
ting eq. (fTU| to reproduce eq. (jHJ precisely @. First, 
\det(M po — I)\ generically does not have the form of a 
simple exponential like p r . Second, no value of \x po in 
eq. (|10p will result in the overall minus sign on the right- 
hand side of eq. ©• 

The generalization of eq. (TTJJ1) to class C has been con- 
sidered by Gnutzmann et al |l4| . As a prototypical class- 
C system, they studied a Fermi sea of electrons (with the 
Fermi surface at E = 0) in a hard- wall billiard in a strong 
magnetic field (to break time- reversal invariance). There 
are then both electron and hole excitations, and 173 is 
defined to be +1 for electrons and —1 for holes. Part of 
the billiard boundary is superconducting, and this leads 
to Andreev reflection: when hitting the superconducting 
boundary, an electron turns into a hole (and vice versa) 
and "retroflects" , initially retracing the incoming path. 
There is an extra phase factor of — i for each Andreev re- 
flection, in addition to the Maslov phase. In general, the 
action of a primitive periodic orbit takes the form [14| 



S po (E) = S$(E) + S$(E), 



(11) 



where S p e J {E) [S P o{E)\ is the action of those segments 
of the orbit where the excitation is an electron [hole] . For 
a given segment, 



SW(E) = -S&J-E). 



(12) 



Gnutzmann et al show that the dominant periodic orbits 
are self-dual. A self-dual orbit includes an odd number 
Na of Andreev reflections, and is traced twice, with each 
segment traced once as an electron and once as a hole. 
For a self-dual orbit, we therefore have 



S po (E) = S$(E) - S$(-E) ~ Et p 



(13) 



where r po — 2dSpJ /dE is the period of the complete 
twice-traced orbit. The Maslov phases of the two tracings 
cancel, but the factor of —i for each Andreev reflection 
results in an extra overall factor of (—i) 2NAT = (— l) r , 
where r is the number of repetitions of the complete orbit. 
Finally, there are two factors of the inverse square-root 
of the stability determinant, one for each single tracing. 
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The final result is therefore [13, EH 



po r 



- r|det(M; o -/)| 



Eq. (fH|) bears a much stronger resemblance to eq. © 
for the Riemann zeros than does eq. (jlOp . The dominant 
orbit actions are linear in E, and the primitive orbits 
contribute with the correct sign. 

We can improve the agreement if we hypothesize that 
the underlying dynamical system has primitive periodic 
orbits that are labelled by both a prime p and another 
integer k = 0, 1, . . . (rather than by a prime p alone), and 
that, for a primitive orbit so labelled, r po — 2 k \ogp and 
|det(M£ - I)| = cxp(rr po /2) 0. With this ansatz, we 
have (setting % = 1) 



AU;(£) = i£E 



(-1)' sin^Tp,,) 



Hdet(M; o -/)| 
1 ( _1 ) r sin(2 fc ri;iogp) 



k=0 r=l 



exp(2 fc r logp/2) 



We now use the mathematical identity [1] 



oo oo 



EE^/(^) = -E> 

fc=0 r=l r=l 

Thus eq. (|T5|) becomes 



(^) = - i EE- 



1 sin(riJlogp) 



n p T^i r exp(rlogp/2) 



(15) 



(16) 



(17) 



which matches eq. ([5]) precisely. Thus, while the even 
repetitions contribute with the wrong sign in eq. (Il~4"l) . 
these contributions can in principle be balanced (in a 
class-C system) by correct-sign contributions from other 
primitive orbits. 

Next we consider our results in comparison with some 
earlier work. 

Connes [HI has suggested that the minus sign in eq. ([5]) 
should be explained by having the Riemann zeros be 
missing eigenvalues in an otherwise continuous spectrum 
of an appropriate hamiltonian H. This would explain 
why all repetitions contribute with the same sign, but 
leaves open the fundamental problem that matching Rie- 
mann zeros to missing eigenvalues does not allow for a po- 
tential proof of the Riemann hypothesis by demonstrat- 
ing that C(^ +iE) oc det(E — H) . Instead, Connes shows 
that the Riemann hypothesis is equivalent to a certain 
trace formula for a hamiltonian with the desired continu- 
ous spectrum. In the present work, we have provided an 
alternative explanation for the sign discrepancy that still 
allows for the original formulation of the Hilbert-Polya 
conjecture. 

Berry and Keating [l9[ have suggested that the quan- 
tum hamiltonian H corresponding to the Riemann zeta 



function should take the form of some quantization, on 
some compactified phase space for one degree of freedom, 
of the classical hamiltonian h(x,p) — xp. Here we note 
that this hamiltonian would be in class D. To see this, 
consider the simplest hermitian quantization on an un- 
compactified phase space, H = \{XP + PX), where X 
and P are the position and momentum operators. If we 
take matrix elements of this hamiltonian between basis 
states with real position-space wave functions, we get a 
hamiltonian matrix of the form H = A, where A is imagi- 
nary and antisymmetric. This characterizes hamiltonians 
in class D 0] . Class- 1? systems have broken time- reversal 
invariance, and hence have eigenvalues with a statistical 
distribution governed by GUE. However, since a class- D 
system does not have the extra nonclassical two-valued 
degree of freedom, eq. (fTU)) for N OBC (E) applies, and so 
the generic sign discrepancy with eq. ([8]) is still present. 

In conclusion, the combination of the Hilbert-Polya 
conjecture (that the imaginary parts of the nontrivial 
zeros of the Riemann zeta function are the eigenvalues 
of some quantum hamiltonian) with the Katz-Sarnak 
conjecture (that the Riemann zeta function is a mem- 
ber of a family of L functions related by a symplectic 
symmetry) implies that a hamiltonian whose eigenval- 
ues are the imaginary parts of the Riemann zeros should 
reside in class C of the Altland-Zirnbauer classification 
scheme. This implies that the hamiltonian should incor- 
porate a nonclassical two- valued degree of freedom. Sys- 
tems in class C generically have primitive periodic orbits 
that contribute to the density of the Riemann zeros with 
the correct sign, further strengthening the argument that 
class C is the right arena to search for the elusive Rie- 
mann hamiltonian. 
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